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Resume 

We report on the experimental measurement of the backscattering matrix of a weakly scattering me- 
dium in optics, composed of a few dispersed gold nanobeads. The DORT method (Decomposition of the 
Time Reversal Operator) is applied to this matrix and we demonstrate selective and efficient focusing on 
individual scatterers, even through an aberrating layer. Moreover, we show that this approach provides the 
decomposition of the scattering pattern of a single nanoparticle. These results open important perspectives 
for optical imaging, characterization and selective excitation of nanoparticles. 

PACS numbers: 42.25.-p,07.60.-j,78.67.-n,42.15.Fr 
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The propagation of waves in inhomogeneous media is a fundamental problem which implies 
applications in many domains of wave physics, ranging from optics or acoustics to solid-state 
physics, seismology, medical imaging, or telecommunications. Conventional focusing and ima- 
ging techniques based on the Born approximation generally fail in aberrating or strongly scattering 
media due to the phase distorsions and the multiple scattering events undergone by the incident 
wavefront. Recent advances in light manipulation techniques have allowed great progresses in 
optical focusing and imaging through complex media. Following pioneering works in ultrasound 
EEL Vellekoop et al. 0) showed how light can be focused spatially through a strongly scattering 
medium by actively shaping the incident wavefront with a spatial light modulator (SLM). Seve- 
ral groups E]-[6) have recently extended this approach to the time domain and demontrated the 
spatio-temporal focusing of femtosecond pulses through highly scattering media. 

A general treatment of such linear complex media is provided by the scattering matrix forma- 
lism, as defined in Q for instance. For ultrasound or microwaves, such a formalism is particularly 
appropriate as the wave field can be generated and detected by arrays of independent elements. Po- 
poff et al. (SI [3 recently extended this method to optics and measured the transmission matrix for 
light between a SLM (emission) and a CCD camera (reception) through a scattering medium. The 
latter is a subpart of the scattering matrix linking the forward scattered field to the incident field. 
The experimental access to this transmission matrix is particularly important since it contains 
many informations on the propagation medium. In particular, it should give access to the open 
channels that allow to transmit light through the scattering sample IfTl fTOlfTTl . 

In this Letter, we study another subpart of the scattering matrix : the backscattering matrix 
K (BM). The knowledge of the BM brings more fundamental insight into the medium. Initially 
developed for ultrasound lfT2l - [T4ll and then extended to microwaves 031 EI), the DORT method 
(French acronym for Decomposition of the Time Reversal Operator) takes advantage of the BM to 
focus iteratively by time reversal processing on each scatterer of a multi-target medium. Derived 
from the theoretical study of the iterative time reversal process, this approach relies on the diago- 
nalization of the time reversal operator KKl A one-to-one association between each eigenstate 
and each scatterer allows focusing on each scatterer by transmitting the computed eigenvectors. 
Technically, a singular value decomposition (SVD) of the BM is used to compute the eigenstate 
of the time reversal operator. In this letter, we will demonstrate the DORT method in the optical 
regime and show its efficiency in weakly scattering media. An experimental setup is proposed 
to record the BM in optics. Using the DORT method, we perform selective focusing on 100 nm 



2 



gold nanoparticles through an aberrating medium, opening an original alternative to conventional 
adaptive optics [17J. In a second part, we show how the DORT method can be used to discrimi- 
nate and selectively excite the different scattering patterns of a single nanoparticle, hence paving 
the way towards a better characterization, imaging and control of nano-objects for instance in 
plasmonics lfT8ll or for super-resolution lfT9ll . 

The experimental setup used for the measurement of the BM is displayed in Fig.[TJ To generate 
the incident wavefront, the beam from a single mode laser at 532 nm is expanded and spatially 
modulated by a liquid crystal SLM. Using a polarizer and an analyzer, we choose an appropriate 
combination of polarizations before and after the SLM to achieve a 27r-phase modulation with 
a residual amplitude modulation below 15%. Light is focused onto the surface of a glass slide 
on which is deposited a dilute solution of 100 nm isotropic gold beads. We image the output 
backscattered field pattern onto a first CCD camera (CCD1) thanks to a beamsplitter. To minimize 
specular reflections, the two polarizers are in a cross-polarized configuration. A reference arm can 
be used to make the backscattered field interefere with a plane wave onto CCD1. The system is set 
such that a pixel of the SLM corresponds approximately to an incident k- vector. An aberrating slab 
is placed 0.5 mm from the image plane between the scattering sample and the objective. It consists 
in an irregular PDMS polymer layer with a mean thickness of 200 fim deposited on a cover glass. 
The CCD1 records the image of the intensity in the plane containing the gold particles, with and 
without the aberrator. For control purposes, another CCD camera (CCD2) images the same plane 
in transmission, through a high NA objective. 

In the experimental setup shown in Fig. [TJ the input and output modes correspond to SLM and 
CCD1 macropixels, respectively. The BM corresponds to a matrix K whose complex coefficients 
k mn connect the optical response (in amplitude and phase) between the n th SLM pixel to the m th 
CCD1 pixel. In other words, the complex optical field on the m th outgoing pixel reads E™* = 
kmnE™ with E™ the field on the n th input pixel. K is of dimension M x N, with M and TV 
the number of pixels contained by CCD1 and the SLM, respectively. The inter-element responses 
k mn are measured following the procedure described in [8], i.e using phase- shifting interferometry 
and a Hadamard input basis. 

A singular value decomposition (SVD) of the BM consists in writing K = USV^ where the 
subscript t denotes the transpose conjugate. A is a diagonal matrix containing the real positive sin- 
gular values Xi in a decreasing order Ai > A 2 > • • • > X N . U and V are unitary matrices whose 
columns correspond to the output and input singular vectors Ui and Vi, respectively. The DORT 
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Figure 1: Schematic of the apparatus. The 532 nm laser is expanded and reflected off a Holoeye LC-R 
2500 liquid cristal SLM. Polarization optics (polarizer Pi and analyser P2) allow phase modulation. The 
modulated beam is focused on the plane containing the gold particules by a 20x objective (NA = 0.5) and 
the backscattered intensity is imaged by a first CCD-camera (CCD1). The second CCD-camera (CCD2) 
images the same plane in transmission with a 60x oil immersion objective (NA = 1.4). Analysers P2 and P3 
are oriented in a cross-polarized configuration. L, lens. BS, beamsplitter. 

method consists in using the SVD of the BM to efficiently detect and separate the responses of 
several scatterers in homogeneous or weakly heterogeneous media lfT3ll . Indeed, under the single 
scattering approximation and for pointlike scatterers, each scatterer is associated mainly to one 
significant eigenstate. The latter is linked to a non-zero singular value that relates to the scat- 
terer reflectivity. In first approximation, each singular vector Vi of K corresponds to the complex 
amplitude mask which, when applied to the SLM, focuses onto the associated scatterer. 

Figs. [2] and [3] show the experimental results of the DORT method applied to the BM matrix 
when several 100 nm gold beads are dispersed on the slide placed behind an aberrating layer. For 
this experiment, we have used the field generated by the unmodulated part of the illumination as 
the reference for the interferometric measurement rather than the plane wave from the reference 
arm. This method detailed in [8] allows for a better stability without being detrimental to focusing 
applications. The numbers of pixels have been fixed to N = M = 1024. Fig. [2]^a) shows the 
image measured by CCD1 in reflection when a plane wave is sent from the SLM. The image of 
the medium is drastically damaged by the aberrating layer : the number of gold beads as well as 
their position and reflectivity cannot be deduced from this image. In contrast, the DORT method 
gives an unambiguous answer to these questions. The singular value distribution (Fig. [2p3)) of 
the BM exhibits two parts. The continuum of low singular values is associated to the noise in 
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FIGURE 2: Experimental results for selective focusing on nanobeads through an aberrating layer using the 
DORT method, a. aberrated image of the object plane on CCD1 in reflection for a plane wave illumination 
b., singular value distribution of the measured BM. c.-f. ( resp. g.-j.) ; images in reflection by camera CCD1 ( 
resp. in transmission by the control camera CCD2), when sending the four first singular vectors. All images 
from CCD1 (resp. CCD2) share the same intensity scale. 



the BM (residual specular reflections, laser fluctuations, CCD readout noise etc.). A few higher 
singular values are associated to the brightest gold nanoparticles since the singular value is directly 
proportional to the reflectivity of the scatterer. The input singular vector Vi corresponds to the 
wave-front focusing on the i th particle and the output singular vector Ui is the corresponding 
scattering pattern measured by CCD1 in reflection (see Figs.[2fc,d,e,f)). Figs.[2jg,h,i,j) correspond 
to the control images measured in transmission by CCD2, when the first four singular vectors 
Vi are back-propagated from the SLM. These images show the selective focusing performed by 
DORT : each singular vector back-propagates selectively on one gold bead without illuminating 
the other ones. Quantitatively, an intensity enhancement of ~ 20 is observed on the four gold beads 
position when DORT focusing is performed compared to plane wave illumination. 

Fig. [3] shows explicitly the effect of aberrations by comparing the singular vectors Vi and Ui 
(associated to the brightest scatterer) with and without the aberrating layer. Figs. [3fb) and (c) dis- 
play the phase of Vi on the SLM in both situations. In free space, the wavefront focusing on 
a particle corresponds to a Fresnel phase zone plates (Fig. ^b)). In contrast, in presence of an 
aberrating layer (Fig. [3jc)), Vi is strongly modified to compensate for the wavefront distortions. 
Figs.|3ja) and (d) display a zoom on the corresponding focal spots recorded by CCD1. The compa- 
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rison of both figures shows the damaging effect of the aberrating layer and highlights the benefits 
of the DORT method in this configuration. 

Figure 3: Phase of the first input singular vector Vi (b, c) and the corresponding focal spots recorded by 
CCD1 (a, d) in free space and in presence of an aberrating layer, respectively. 

So far, we made the assumption that each scatterer is only associated to one single singular 
value. However, this is not strictly true : there can be several singular values associated to a single 
scatterer, each one corresponding to a different scattering pattern of the particle Il20l - [22l . To in- 
vestigate the different eigenstates associated to one nanoparticle, we now study analytically and 
experimentally the DORT analysis of the BM without any aberrating layer and with only one gold 
particle in the observation area. 

Let us assume that the incident field is linearly polarized along x, z is the propagation axis and 
y the polarization observed in reflection (cross-polarized configuration). The BM can be decom- 
posed into the product of three matrices : K = K".a.K'. Here K' is a 3 x TV matrix describing the 
propagation of light from the TV pixels of the SLM to the scatterer. The three lines of K' corres- 
pond to the different orientations of the electric field at the particle position, a is the polarisability 
tensor of the particle linking the local electric field to the dipole moment p of the particle. K" 
is a M x 3 matrix which connects the dipole moment p to the y-component of the electric field 
measured on the M pixels of CCD1. We will assume that the particle is isotropic : a — al with a 
the particle polarizability and I the identity matrix. 

The electric field illuminating the particle in the focal plane can be deduced from [[23l [24] and 
leads to the following expression for the coefficients of K' : 

k' l3 = / o (ej,p fl )+/ 2 (ej,p fl )cos(20 fl ) (1) 

s v ' 

negligible 

k' 2j = / 2 (0],p s )sin(2^) (2) 
k' 3j = 2j/ 1 (flj,p s )cos(0 s ) (3) 

4> s and p s are the polar coordinates of the particle in the focal plane (see Fig. |4^a)). 6>] < Q™ ax is 
the angle made by the k-vector associated to the j th element of the SLM with the optical axis. The 
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functions J , h and 7 2 are given by : 

J o (0,p) = Vcos 9 sin 0(1 + cos 6) J (2np sin 6) (4) 

7i(0,p) = V / cos0sin 2 0J 1 (27rpsin0) (5) 

h(0,p) = Vcos#sin#(l -cos#)J 2 (27rpsin#) (6) 



with «/i and J 2 the zero, first and second order spherical Bessel functions of the first kind, 
respectively. 

The ^-component of the backscattered field can be also deduced from the calculations of ll23ll 
and leads to the following expression for K" : 

Qrnax 

k' k \ = / d6 2 I 2 (6 2 ,p k )sm(2<f> k ) (7) 
Jo 

Qmax 

K2 = I dd 2 [I (e 2 ,Pk) + 12(62, Pk)coa(2<l> k )] (8) 
Jo v v ' 

negligible 

Qmax 

Ks = 2 j 2 d9 2 h(9 2 ,p k )sm((f> k ) (9) 
Jo 

Pfe and (j> k are the polar coordinates of the k th pixel of CCD1 and 6 2 is the angle made by a beam 
diffracted by the particle with the optical axis before going through the microscope objective (see 
Fig.[4|b)). Hence, the BM depends on the following parameters : the position of the particle in the 
focal plane (0 s ,p s ), the numbers of SLM and CCD elements (N and M) and the input and output 
maximum angles (0™ ax and 6™ ax ). 6™ ax is directly given by the numerical aperture of the optical 
system, but Q™ ax can be smaller if the area of illumination is smaller than the area of the pupil of 
the microscope objective. Theory predicts that the x 9 y and z components of the particle dipole 
moment p give rise to different scattering patterns : monopolar for y (Eq{8]), bipolar for z (Eqj9j) 
and quadrupolar for x (EqJTJ) (see Fig|4|d)). Those scattering patterns are strictly orthogonal if 
they are observed over a solid angle Q = 2tt. Hence, in theory, the output singular vectors Ui of K 
should correspond to the different scattering patterns of the nanoparticle (i.e the columns of K"). 
In practice, this is not strictly true since only a limited part of the half space is accessible in our 
experiment. 

Our theoretical predictions are now confronted to the experimental results. As an exact measu- 
rement of the output singular vectors is needed, the complex wave field has to be measured directly. 
A reference arm is now used to make the backscattered wave interfere with a reference plane wave 
on CCD1 (see Fig{T]). The output singular vectors Ui correspond then to the complex scattering 
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FIGURE 4: Singular value decomposition of the experimental and theoretical BM with a single scatterer. (a) 
and (b) represent the emission and the recording geometry of the experiment. (c) and (d) show the singular 
value distributions of the measured and calculated BM. The three first output singular vectors are shown in 
the insets. 

patterns of the particle. The number TV of controlled segments on the SLM has been lowered to 16 
for stability contraints over the time of the experiment. The result of the DORT analysis is shown 
in Figj4| The experimental and theoretical singular values distribution as well as the three first out- 
put singular vectors of the BM are compared in Figj4jc) and (d), respectively. The experimental 
values of Q™ ax and 6™ ax have been used for our theoretical calculations (0™ ax = 22°, 6™ ax = 68° 
). The position of the particle in the focal plane has been arbitrarily fixed (cj) s = 7/87T, p s = 665 
nm) such that the theoretical and experimental ratios between the first and the second singular 
values are equal.Our theoretical calculations predict that only the three first singular values are 
non-zero, each one corresponding to the independent scattering patterns of the three orthogonal 
dipolar radiations (along x, z and y resp.). The experimental results display a continuum of low 
singular values corresponding to noise and three higher singular values (Fig|4jd)). The two first 
output singular vectors match the scattering patterns of the x— and z— components of the nano- 
particle dipole moment, as predicted theoretically. The third output singular vector has a geometry 
comparable to the second one, but is more affected by the experimental noise since a significant 
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part of its energy is observed outside the particle. This singular vector might correspond to a cou- 
pling between the second singular state (z— component) and eigenmodes of the noise background. 
The scattering pattern associated to the y— component of the dipole moment is drowned in noise 
and could not be observed. 

The DORT method allows to excite independently each component of the nanoparticle dipole 
moment : such an experiment would require to increase the number of SLM segments for quan- 
titative results and a near field scanning optical microscope (NSOM) to probe the near-field of 
the nanoparticle. For larger scatterers, the singular spectrum may be even richer : more radiation 
modes would be excited and a better characterization of the nanoparticle could be obtained (size, 
position, permittivity etc.). In particular, one could discriminate the quadrupolar radiation of the 
particle, which may open important perspectives in plasmonics. 

In summary, we have performed an experimental measurement of the BM of a weakly scatte- 
ring medium. We have used this BM to demonstrate the DORT method in optics and proved that 
it allows selective focusing on nanoparticles and their caracterization. This approach could lead to 
many applications in nano-optics, sensing or imaging and entails the interest of the BM in optics. 
We are confident that measuring the BM in more complex media will raise interesting problems in 
optics. For instance, the measurement of the BM in highly scattering media may allow the sepa- 
ration of the single and multiple scattering contributions ll25ll as well as the study of the coherent 
backscattering cone and of the diffusive halo |[26l . 
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